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We analyze several product measures in the space of mixed quantum states. In particular we 
study measures induced by the operation of partial tracing. The natural, rotationally invariant 
measure on the set of all pure states of a N x K composite system, induces a unique measure in the 
space of N X N mixed states (or in the space of K x K mixed states, if the reduction takes place with 
respect to the first subsystem). For K — N the induced measure is equal to the Hilbert-Schmidt 
measure, which is shown to coincide with the measure induced by singular values of non-Hermitian 
random Gaussian matrices pertaining to the Ginibre ensemble. We compute several averages with 
respect to this measure and show that the mean entanglement of N x N pure states behaves as 
<N ■ InJV- 1/2. 
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I. INTRODUCTION 



Recent developments in the emerging field of quantum information increased an interest in investigation of the 
properties of the set of density matrices of a finite size. To characterize quantitatively properties of typical density 
matrices it is necessary to define a certain measure /z in this set, which will enable one to compute the desired 
averages over all N x N density matrices with respect to this measure. Only having chosen the measure it makes 
| sense to ask about the probability to find a random density matrix representing an entangled state or to compute 
the average entanglement of mixed quantum states (lp|| . Alternatively one may search for a statistical ensemble of 
■y | density matrices which corresponds to minimal prior knowledge about a quantum system [QJ^]. 

In the space of pure states of an N dimensional Hilbert space H.n, isomorphic with the complex projective space 
CP W_1 , there exist a unique 'natural' measure, Pn(\4>)), induced by the Haar measure on the unitary group U(N). 
q-i In other words, the random pure state, defined by action of a random unitary matrix on an arbitrary reference state, 
|0) = U\<fio), may be represented (in an arbitrary basis) as a certain row (or column) of the random matrix U. The 
. J_h | measure Pjv(|0)) = Pn{U\<P)) in the space of pure states is thus distinguished by the rotational invariance of the Haar 
^ ■ measure [||-||. 

This property of the set of pure states may suggest that there exists a natural measure in the space of mixed 
quantum states Mn, distinguished in an analogous way. However, this might not be the case: it is not at all simple 
to select a single measure in M. n and provide rational arguments for its superiority • Any density matrix g 

is Hermitian and may be diagonalised by an unitary matrix U. It is thus tempting to assume that the distributions 
of eigenvalues and eigenvectors of g are independent, so fi factorizes into a product measure A x v. For the measure 
v on the unitary group, which determines the statistics of the eigenvectors forming U, one may take the unique Haar 
measure on U (N) . On the other hand, it is much more difficult to find convincing arguments, which allow one to pick 
out the unique measure A defined on the simplex of eigenvalues. 

To consider a possibility of defining induced measures in the (N 2 — 1) dimensional space Al at it is worth to discuss 
the procedure of purification. Any mixed state g, acting on TLn, may be purified by finding a pure state |\&) in the 
composite Hilbert space Hn ®1~Ini such that g is given by the partial tracing over the auxiliary subsystem, 

|*) -^£? = Tr 2 (|*)(*f) . (1.1) 
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In a loose sense the purification corresponds to treating any density matrix of size N as a vector of size N 2 . For a 
given g finding a corresponding pure state | "J/) does not have a unique solution since 

= Tr 2 (|*)(tf|) =Tr 2 (V)(*'l) where |*') = (I ® I7 2 )|tf) (1.2) 

where C/2 is an arbitrary unitary matrix of size N acting locally on the auxiliary subsystem. 

It is now simple to formulate the idea of the induced measure: we choose the natural measure P jy2 ( 1^)) in the space 
of purified states, and look for the measure Pn,n(q) induced by the operation of partial tracing (1.1). This idea, put 
forward by Braunstein Q , was further developed by Hall Q . He considered the ensemble of pure states in the NK 
dimensional space TIn <8> 'H-k distributed according to the the natural measure Pnk(\^))- The partial tracing over 
the variables of TIk defines uniquely the induced measure Pn,k{q) in the space of mixed states Mn- Hall found 
explicitly the probability distributions P2,k{b) and computed mean entropy averaged over these measures. The aim 
of this work is to extend his results providing the general formula for Pn,k(q)- Moreover, we propose other methods 
of defining induced measures via the Hilbert-Schmidt space and establish a link between them. 

The space of linear operators acting ouHn, equipped with the scalar product (A\B) =tvA^B, is called the Hilbert- 
Schmidt space, Tins- It has N 2 complex dimensions and may be represented by the algebra of complex matrices of 
size N. Any non-zero matrix A may be projected into the space of density matrices by |lO[ 

A /4t 

It is easy to see that AA^ is Hermitian and positive (strictly speaking, non- negative), while the rescaling by the 
norm assures the trace condition trp=l. Obviously one may consider an alternative, entirely equivalent projection 



A — > g' = A^ A/trA^ A, which agrees with (1.3) for normal A. 

Two matrices A and A 1 produce the same state g, if there exist an unitary V, such that A' — AV. An arbitrary 
measure P(A) in the Hilbert-Schmidt space induces thus a certain measure in the space M. of density matrices. The 
eigenvalues of g are equal to the rescaled squared singular values of matrix A, where SV(A) = ^/cig(AA^) [|l). Note 
that the singular values of A and A' are the same. For P(A) one may take the measure of Ginibre ensemble of 
(non-Hcrmitian) random Gaussian matrices. In this paper we prove that the measure induced in this way is equal to 
Pn,n{q), which is also shown to coincide with the measure related to the Hilbert-Schmidt metric. 

Our paper is organized as follows. In Section II we write down the natural measure in the space of pure states and 
provide a review of different measures used in the space of the density matrices. Section III is devoted to measures 
induced by partial tracing. In section IV we discuss possibilities of generalizing our results and provide certain families 
of measures in M.. The paper is concluded in section V. 

II. MEASURES IN THE SPACE OF DENSITY MATRICES 

A. Natural measure on the space of pure states 

Before concentrating on the set A4n of the density matrices of size N, let us discuss the space of pure states. For 
concreteness we start our considerations with the four dimensional Hilbcrt space - the simplest case important from 
the point of view of quantum entanglement. The set of pure states of a 4 dimensional Hilbert space forms a complex 
projective space CP 3 , on which the natural uniform measure exists. To generate random pure states according to 
such a measure on this 6 dimensional manifold, we take a vector of a random unitary matrix distributed according to 
the Haar measure on t/(4). The Hurwitz parametrisation ||l2]| gives 

I*) = (cos^^in^cos^e^^in^sin^cos^ie^^in^sin^sin^ie^ 1 ), (2- 1 ) 

where $ k G [0, tt/2], and <p h G [0, 2tt) for k = 1, 2, 3. 

A uniform distribution over almost all of CP 3 is obtained by choosing the uniform distribution of the 'polar' angles; 
P{fk) = 1/2-7T. In analogy to the volume element on the sphere, the 'azimuthal' angles § k should be taken in a 
nonuniform way, with the probability density p"2] ] 

P{d k ) = ksm(2d k )(smtf k ) 2k - 2 for $ k G [0, tt/2], k= 1,2,3. (2.2) 
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In practice it is convenient to use auxiliary independent random variables £fc distributed uniformly in [0, 1] and to 
set $fc = arcsin(^^ 2fe ) . In order to generalize the parametrisation (2.1) for a vector belonging to the TV-dimensional 



Hilbert space, one uses (N — 1) pola r an gles uniformly distributed in [0, 2tt), and (N — 1) independent azimuthal 
angles distributed according to (|2.2|) |13|. The volume element in space CP N ~ X reads thus 



N-l 



dv = Y[ costfkismtfkf^dtfkdipk, (2.3) 



fc=i 

while the total volume of the manifold of pure states 

N-l f7 ,/2 N-l r 2n 



Vn=U cos^^sin^) 2 "- 1 ^ 11 / dy> fc = fjv _ lV ( 2 - 4 ) 

fc=i ^° fe=i ^° ' 



is proportional to the product of the volume of the (N — l)-d torus, (built of the phases (fk), times the volume of 
the (N — l)-d simplex (built of the angles ?9fc) |H. The natural measure (2.3) in CP^ -1 induces the uniform 



measure in the simplex Sn ]3j . In appendix A we show that the same measure is given by squared absolute values of 
N independent random complex Gaussian numbers Zj, rescaled as |zj| 2 — > \zi\ 2 / i \ z i\ 2 ■ Thus a random pure state 
may be constructed of N appropriately rescaled complex random numbers drawn according to the normal distribution. 

It is worth to note that the natural measure ( [2.3| ) is compatible with the Fubini-Study metric on CP^ -1 , in the 
sense that the measure of a certain e- neighbourhood of any chosen point is the same j[o| . The average value of some 
function of a state, X^ = X(\ijj)) (e.g. an observable (X)^,), over the entire manifold of pure states is then defined as 



X = -i- / X^dv{\i))). (2.5) 
'" /cp"- 1 



V N 



B. Product measures in the space of mixed states 



Consider a Hermitian, nonnegative matrix g of size N, normalized as trg = 1. It represents a mixed state belonging 
to the space Mn and may be diagonalized by a unitary rotation. Let B be a diagonal unitary matrix. Since 

g = UAU^ = UB\B*U\ (2.6) 

the rotation matrix U is determined up to N arbitrary phases entering B. The number of irrelevant phases increases, 
if the spectrum of g is degenerated (see e.g. p5|). The total number of independent variables used to parametrize a 
density matrix g is equal to N 2 — 1, provided, no degeneracy occurs. In recent papers discussing the relative volume 
of the set of entangled states [0J|] , we considered the product measures 

/i = A x v. (2.7) 

The first factor defines a measure in the (N — 1) dimensional simplex of eigenvalues entering A, since due to the trace 
condition 52»=i ^ = 1- The measure u, defined in the space of unitary matrices U(N), is responsible for the choice of 
eigenvectors of g. It is natural to take for v the Haar measure v u on U(N), so that the probability distributions are 
rotationally invariant, P(g) = P(U gU^). On the other hand, as pointed out by Slater Q, there exist several possible 
choices for the measure A. 



C. Measures induced by random vectors 



In Rcf. ||] we defined a unitary product measure by taking for the vector A = {A^} the squared moduli of complex 
elements of a column (say, the first column) of an auxiliary random unitary matrix V drawn with respect to the Haar 
measure v u on U (N) 

A 4 = |^i| 2 . (2.8) 

The distribution P U (A) obtained in this way is nothing else but the distribution of components of a random vector 
on CP N ~ 1 , uniform on the simplex of eigenvalues Q|. In a similar way one defines the orthogonal product measure 
where Xi = |Oa| 2 and O is an orthogonal matrix drawn with respect to the Haar measure v on O(N). 
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Let us stress that the name of the product measure (orthogonal or unitary) is related to the distribution A on the 
simplex of eigenvalues, while the random rotations U are always assumed to be distributed according to the Haar 
measure v u in U(N). Thus both measures may be linked directly to the well known Gaussian unitary (orthogonal) 
ensembles of random matrices Cq,0|, referred to as GUE and GOE. The measure /i Q is determined by squared 
components of a real eigenvector of a GOE matrix, while the measure fi u may defined by components of a complex 
eigenvector of GUE matrices. 

It is easy to see, that these measures on the simplex are special cases of the Dirichlet distribution [|| 

A^Ai.-.-.A^v-i) =a s Ar 1 -..A^ 1 1 (l-A 1 -. ..- \n-iY~ 1 - (2.9) 

where s > is a free parameter and a s stands for a normalization constant. The last component is determined by 
the trace condition Aat = 1 — Aj. — • ■ ■ — Ajv-i' The uniform measure Ai describes the unitary measure fi u = Ai x u u , 
while the case s = 1/2 corresponds to the orthogonal measure, ii = A±/2 X u u ||. The latter is related to the Fisher 
information metric [Q , the Mahalonobis distance pO] and Jeffrey s' pr ior |2~T 



It is instructive to consider the limiting cases of the distribution (2.9). For s — > one obtains a singular distribution 
peaked at the pure states [0]) while in the opposite limit s — ► 00, the distribution is concentrated on the maximally 
mixed state described by the vector A = {1/N, . . . , 1/N}. Changing the continuous parameter s one can thus 
control the average purity of the generated mixed states. 



D. Measures related to metrics 

There exist several different distances in the space of mixed quantum states - (see e.g. |2^,|23|]). Each metric 
generates a measure, since one can assume that drawing random matrices from each ball of a fixed radius is equally 
likely. The balls are understood with respect to a given metric. We shall consider here two most important examples: 
the Hilbert-Schmidt metric and the Bures metric. 

The Hilbert-Schmidt distance between any two density operators is given by the Frobenius (Hilbert-Schmidt) norm 
of their difference 



Dhs(qu Q2) = VWoi - 92) 2 }- (2.10) 

With respect to this metric the set of the two levels mixed states has the geometry of the Bloch sphere and its interior. 
For any dimension N the infinitesimal distance takes a particularly simple form 

(ds HS ) 2 = Tr[(Sg) 2 } (2.11) 

and allows one to compute the volume element, normalize it by appropriate constants Cm and to obtain the probability 
distribution 

N N 

Phs(Xu ■ • • , Ajv) = C N 6(1 - J2 x i) Y[( x i ~ x k) 2 - (2.12) 

i—l j<k 

This joint probability distribution, derive d by Hall Q, defines the Hilbert-Schmidt measure &hs m the space of 



diagonal matrices, and a product measure (2.7) in the space of density matrices 



Another measure is related with the Bures distance in the space of mixed quantum states, |24 



Db(qi,Q2) = V2(l-tr[(^ /2 2 f?i /2 ) 1/2 ]), (2-13) 



It corresponds to the minimal Fubini-Study distance between all possible purifications (1.1) of both mixed states Qi 
and Q2 ■ With respect to this distance the 3-dimensional set of N = 2 density matrices exhibits the geometry of a half 
of the hypersphere S 3 € R 4 , with pure states at the hyper- meridian p6| ]. 

The infinitesimal Bures distance between g and Sg was found by Hiibner p7| 



j,k=i 
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where Xk and |fc) represent eigenvalues and eigenvectors of g. Hall computed the corresponding volume element, and 
received the Bures probability distribution M 

P B (X U ...,X N ) = C' N f-^$ n ( 2 - 15 ) 

This measure was later analyzed by Slater [T^ ], who found the normalization constants Cy for N — 3,4 and 5 and 
together with Byrd analyzed the geometry induced by the Bures measure for N = 3 |2q |. The Bures measure is 
monotone, i.e. it does not increase under the action of completely positive, trace preserving maps p2fl. 



E. TV = 2: Measures in the Bloch ball 



To experience some features of different measures consider the set of 2 x 2 density matrices. Such a matrix can be 
expressed by the Pauli matrices a as g = 1/2 + a • r, where r is a three dimensional vector of modulus r < 1/2. Pure 
states correspond to the Bloch sphere, r = 1/2, while the center of the ball, r — 0, corresponds to the maximally 
mixed state proportional to the identity matrix. The eigenvalues of g read Ai = 1/2 + r, and A2 = 1/2 — r. The 
rotationally invariant Haar measure v u implies in this case the uniform distribution of the pure states on the Bloch 
sphere. The product measures may differ by a different radial distribution P(r). 

The unitary measure Ai is uniform in radius, Pu(r) — 2 for r G [0, 1/2], while the orthogonal measure A]y 2 defined 



by ( [2.9] ) gives the "cosine distribution", P {r) = A/(tt\J\ — Ar 2 ). The Hilbcrt-Schmidt measure ( [2.12 ) reads for N = 2 

P HS (r) = 24r 2 for re [0,1/2] (2.16) 

The quadratic factor relates to the Jacobian in 3-D spherical coordinates, so the Hilbert- Schmidt measure corresponds 
to the uniform coverage of the entire Bloch ball . 

On the other hand, the Bures measure (2.15) implies the following radial distribution, obtained by substituting 
Ai = 1/2 + r and A 2 = 1/2 - r 



P B {r) 



32r 2 



71VI - 4r 2 



(2.17) 



where r G [0,1/2]. As shown in Fig. la this distribution is more concentrated on high purity states with large r, 
located close to the Bloch sphere. 

To analyze the simplest measure P 2 . 2 (g) induced by partial tracing, consider pure state = (x,y,z,t) of the 
2x2 = 4 dimensional Hilbert space H.4. The corresponding density matrix reads 



l*X*l 



xx xy 

* * 

yx yy 

zx zy 

tx* ty* 



yz* yt* 
zz* zt* 
tz* tt* 



Performing the partial tracing with respect to the second subsystem one gets the reduced density matrix 



g A = Tr B {g) 



xx* + yy* xz* - 


Yyt* ' 




b w 


zx* + ty* zz* - 


Ytt* 




w* 1 — 6 



(2.18) 



where b = xx*+yy* and w — xz*+yt* . Its eigenvalues might be written as Ai j2 = \^- r , where r = \ \J{^ — 26) 2 + A\w\ 2 
It has been alread y pro ved by Hall that the induced distribution P 2 . 2 (r) = 24r 2 and coincides with the Hilbert- 
Schmidt measure ( |2.16 ). In the next section we demonstrate that this relation holds for arbitrary N. 

The pure states l^) e H4 may be represented by the Schmidt decomposition |^) = cosa|00) +sina|ll), where 
{|00), 1 1} , 1 10} , |11)} denotes the standard basis in H4 a nd a G [0, 7r/4]. The corresponding eigenvalues of the reduced 
density matrix are Ai = cos 2 a and A 2 = sin 2 a, so Eq. ( [2.16] ) implies the following distribution for the Schmidt angle 



P(a) = 3cos(2a) sin(4a). 



(2.19) 
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For 2x2 systems (N = 2) all measures of pure states entanglement are equivalent and may be expressed as a function 
of the Schmidt angle. For example, the t angle reads 0] r = AdetpA = 4A1A2 = sin 2 (2a), while concurence C is equal 



to yfr = sin(2a). Changing variables in (2.19) we obtain the distributions 



P(t) = -^fi^T and P{C) = ZCy/l - C 2 . (2.20) 
consistent with numerical results of Kendon et al. 63). First moments of these distributions give the mean tangle 



(t) CP 3 = 2/5 and the mean concurrence (C)cp3 = 37r/16 ~ 0.589049, averaged in the sense of (2J5) over the natural 
measure on the space of TV = 4 pure states. 

In a similar way it is straightforward to compute the mean entanglement E. It is just the mean von Neumann 
entropy, H(q) — —trying, of the reduced density matrices 

/•1/2 1 
(E) CP3 = (E)hs = [-rlnr-(l-r)ln(l-r)]24r 2 dr= -. (2.21) 
Jo 

In the natural units we use the maximal entanglement E ma]t — In 2. If one expresses the entanglement E in units of 
In 2 then E max /\n2 = 1 and the mean value reads (E)/\n2 = 1/(3 In 2) w 0.480898. For other orthogonal, unitary 
and Bures measures we obtain, respectively, (E) = 21n2 — 1 w 0.386, (E) u = 1/2 and (E)b = 21n2 — 7/6 w 0.219. 
In the same way we compute the mean inverse participation ratio R = [trg 2 ] -1 of the reduced matrices. This quantity 
may vary from R = 1 (for pure states) to R — N (for the maximally mixed state), and has an interpretation as 
the effective number of pure states in the mixture. For density matrices distributed according to H-S measure one 
has Rhs — 5/4. Performing an analogous integral for the Bures measure we obtain Rb = 8/7, which reveals high 
concentration of this measure on the states of higher purity - see Fig. la. On the other hand, for orthogonal and 
unitary measures we find R Q = 4/3 and R u — 3/2. 




FIG. 1. Distribution of the eigenvalues of density matrices of size N = 2: (a) unitary product measure (uniform distribution 
- dotted line), orthogonal product measure ("cosine distribution" plotted with dashed line), Hilbert-Schmidt measure (solid 
line parabola) and Bures measure (dash-dotted line); (b) measures P-2,,k induced by partial tracing; solid line represents H-S 
measure i.e. K = 2, dashed line K = 3, dash-dotted line K = 4, dotted line K = 5. 



III. MEASURES INDUCED BY PARTIAL TRACING OF COMPOSITE SYSTEMS 



A. Joint probability distribution 



Consider a bipartite N x K composite quantum system. Pure states of this system |\E') may be represented by 
a vector = {^i, ■ ■ ■ , ^m} t , where M — NK is the dimension of the composite Hilbert space Ti = Ha <8> Hb- 
In view of the operation of partial tracing it is convenient to work in a product basis \i) ® \k), where \i) € Ha', i = 
1,...,N while \k) € 7is;fc — 1, . . . , K. A pure state |\E') is then represented by a rectangular complex matrix 
tpih. The normalization condition, ||^|| 2 = tri/jip'' = 1, is the only constraint imposed on this matrix. The density 
matrix acting on the composite Hilbert space Ti., is represented in this basis by a matrix labeled by four indices, 
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Qik,i'k' = (*| ® {^\qW) ® N') = ''Pikipi'k'* ■ The partial tracing with respect to the K dimensional subspace TLb gives a 
reduced density matrix of size N 



K 

g A = Tr B g where g A = ^ ^iktpjk* , (3.1) 

k=l 

while a symmetric operation of partial tracing over the first subsystem leads to a reduced density matrix g B = TiaQ 
of size K, where g B t — ^2%=i ^ik^tt* ■ The natural measure in the space of M = KN dimensional pure states induces 



then the measure Pn,k{q) in the space of the reduced density matrices (3.1) of size TV. Note that the problem 
of measures induced by partial tracing may be considered as a projection of the (KN — 1) dimensional simplex of 
eigenvalues Skn-i into simplices of smaller dimensions 5jv-i or Sk-i- 

Without loss of generality we assume K > N, then g A is generically positive definite. In the opposite case K < N 
the reduced density matrix g A has N — K zero eigenvalues, but the reduced matrix g B of size K is positive definite and 
has the same positive eigenvalues. In any case we are interested only in the distribution of the positive eigenvalues. 
Let us call the corresponding positive reduced density matrix again g — tpip' where ip is now considered as a N x K 
matrix. First we calculate the distribution of matrix elements 



P(g) oc / [dip] S(g - ipip^) <5(Tr<^ f - 1) (3.2) 

where the first delta function is a delta function of a Hermitian Matrix and in the second delta function Tiipip' may 
be substituted by Trg. Since g is positive definite we can make a transformation 

i/j = y/ffip, [#] = detg K [dip] . (3.3) 
Note that [dtp] includes (alternating) factors dipik and dipik* ■ The matrix delta function may now be written as 

5(^(1 - ffi)V0) = A ^9~ N S(l - ffi) , (3.4) 

where the first factor on the right hand side is the inverse Jacobian of the corresponding transformation. As the result 
the distribution of matrix elements is given by 



P(g) oc 9(g) 6(Tr<? - 1) detg 



K-N 



(3.5) 



where the theta function assures that g is positive definite. It is then easy to show by the methods of random matrix 
theory that the joint density of eigenvalues A = {Ai, A2, . . . , Ajv} of g is given by 



Pn,k(X u ...,X n ) = C n , k 5(1 - J2 A *) II X f~ N °^) II( A < - A ^ 



(3.6) 



The square of the Vandermonde determinant results from integrating out the eigenvectors and the normalization 
constant ( see also (4.7)) may be expressed in terms of the Euler Gamma function T(x) 



Cn,k = 



T(KN) 



n^" lr (^-i)r(AT-i + i) 



(3.7) 



written here for K > N . Observe that for K — N the measure Pn,n induced by partial tracing of the pure states in 
(£pjv 2 -i confides with the Hilbert-Schmidt measure ( 2.12 ), so Cn.n — Cjv- Moreover, it equals to the probability 
distribution of squared singular values af of matrices pertaining to the normalized Ginibre ensemble of complex 
matrices A, where ai are given by the square roots of non-negative eigenvalues of AA^ . This is even valid if A is a 
N x K matrix (K > N). To clarify this point instead of starting from equation (3.2) we could have been started from 
a distribution 



P(g) oc / [dip] 5(g - ipipl /Tripipi) exp(-Tr^ T ) 



(3.8) 



Introducing a delta funct ion for t = Tripip* 1 and rescaling ip with \/t we see that both distributions ( |3.2| ) and (3.8) are 
equivalent and given by (3.5). This follows also from the fact that N 2 complex elements of the random matrix A of 
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the Ginibre ensemble may be treated as a random vector distributed according to the natural measure on CP N _1 - 
see Appendix A. 

In this way one may give a si mple recipe to generate a random density matrix according to the HS measure and 
similarly the generalisation (3.5): 

a) prepare a square complex random matrix A of size N pertaining to the Ginibre ensemble (with real and imaginary 
parts of each element being independent normal random variables) , 

b) compute the matrix H = AA* /(trvL4J), (generically positive definite). 

In principle one could proceed further diagonalizing H and writing the density matrix as g — UdW where d is the 
vector of eigenvalues of H and U is a ra ndom unitary rotation drawn according to the Haar measure v, so g is taken 
according to the product measure (2.7). However this step is not necessary, since the Ginibre ensemble is invariant 
with respect to unitary rotations, so one may set g = H. 

The explicit form of the distribution (3.6) for N = 2 was already derived by Hall One could speculate, the 
larger dimension K of the auxiliary space TLb, the more uniform is the induced distribution Pjv,k"(A). This is not 
the case [Q, and the larger K the more is the induced distribution concentrated in the center of the Bloch ball - 
see Fig. 2b. A similar effect is shown in Fig. 2 for N = 3. It shows the probability distributions in the eigenvalues 
simplex - the equilateral triangle situated in the plane A3 = 1 — Ai — A2. The corners represent three orthogonal pure 
states, while the edges denote the matrices of rank 2. The more mixed the sta te, t he closer it is to the center of the 
triangle, which represents the state g* = 1/3. Unitary product measure, Eq. ( |2.9| ) with s = 1, covers uniformly the 
entire triangle, while the orthogonal product measure, s = 1/2, distinguishes states of higher purity - see Fig. 2a. The 
Bures measure ( |2.15 ) with normalizatio n con stant C' N = 35/7T is even more localized at states close to pure - see 
Fig. 2b. The Hilbert-Schmidt measure ( [2.12 ), equal to P 33 (A), is shown in Fig. 2c, while other induced distributions 
(3.6) are plotted in Fig. 2d-2f. Observe that the larger dimension K, the more the distribution P 3 #-(A) shifts in the 



eigenvalues simplex toward its center. Due to the factor (Ai — Xj) in (3.6) the degeneracies in spectrum are avoided 



which is reflected by a low probability (white colour) along all three bisetrices of the triangle. On the other hand, the 
distribution P 3 ,2(A) is singular, and located at the sides of the triangle, which represents the density matrices of rank 
2. It is equal to P2 j3 (A) represented by dashed line in Fig lb. 
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FIG. 2. Distribution of the eigenvalues of density matrices of size N = 3 in the simplex of eigenvalues (a) orthogonal 
product measure, (b) Bures measure, (c) Hilbert-Schmidt measure equal to ^3,3, other measures induced by partial tracing: 
(d) P 3 ,4, (e) P 3 ,5 and (f) P 3>6 . 



B. Averages over induced measure 



Let us now calculate some moments with respect to the Hilbert-Schmidt measure Pj/s(A) = P7v,tv(A). Taking 
Fourier representation of the delta function 6(1 — ^ i Ai) one sees that averages which are homogeneous functions of 
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some Afe are related to the corresponding moments of the Laguerre ensemble: 

P L (X U ...X N ) oc exp(- \) II A f ~ N0 ( x i) li^ - X of- (3-9) 



Here simple moments can be calculated by the method of orthogonal polynomials |L6|. For example one finds for 
K = N and v > -1 

,00 N-l 

(Tr Q v ) H s = B N , V \ dxx v J2 L n (xfe- X , (3.10) 



with the Laguerre polynomials defined by 



L n {x)=e x -(^) n x n e- x (3.11) 
n! ax 



and Bn,v a constant which results from the map of the Hilbert Schmidt measure to the Laguerre measure: 

B N ,u = T(N 2 )/T{N 2 + (3.12) 

Note that we also obtain an explicit expression for the entropy H(g) = — (Trying) by taking the derivative of Eq.( |3.1C| ) 
at v — 1. This can be written as a sum of Digamma functions. For large N it displays asymptotic behaviour 

(H(Q)) HS ^hxN-l+0(^.). (3.13) 

This formula gives one of the main results of this paper - the mean entanglement of pure states of the N x N bipartite 
system, averaged over the natural measure on the space of pure states. It is equivalent to (H(q))hs ~ ln(0.6065iV), 
which stays in agreement with recent numerical results of Lakshminarayan pl[ , who computed mean entanglement 
of pure states of an exemplary chaotic quantum system - the four dimensional, generalized standard map without 
time reversal symmetry. Our result is somewhat similar to the asymptotics of the mean entropy of complex random 
vectors (i3 r (|</)))) CP iv-i sa lniV - 1 +7 @], where the Euler constant 7 - 0.5772. 
The recurrence relation 

xL n {x) = -(n + l)L n+1 {x) + (1 + 2n)L n (x) - ni„_i(a;) (3.14) 

implies for v = 2 

N-l 

{Trg2)HS = ivWTT) ^ {{n + 1)2 + (1 + 2,1)2 + " 2) (3 - 15) 

Finally the sum over n can be done with the result: 

IN 

(Trg 2 ) HS = (3.16) 

Using the recurrence relation for the L n (x) one can find the term oc L n (x) in the Laguerre expansion of x l 1 L n (x) for 
large n. From this we find the following asymptotic behavior for large N: 



This implies for the entropy an asymptotic behavior ( ^.13[ ). Inspecting the recurrence relation further one sees that 
the moments in question are always given by the ratio of two polynomials of N. With this knowledge it is easy to 
calculate some further moments: 
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and 



4 14iV 3 + 107V 

{TTg )HS - (iV2 + l)(JV2 + 2)(iV2 + 3)- (3 - 19) 

Finally let us mention that the result for the second moment can be generalised for density matrices distributed 



according to (3.6) for arbitrary N and K 



The calculation leads to the generalised Laguerre polynomials. Due to the symmetry in K and N it is even not 
necessary to mention what is the largest of both dimensions K or N. The participation ratio reads thus R = 
(NK + 1)/(N + K) and is consistent with recent results of Zanardi et al. |32). They computed the linear entropy 
1 — trp 2 for an ensemble of mixed states obtained by partial tracing of pure states U(\ipi) <g> \ip2)), where \ip{) and 
are pure states distributed according to the natural measure on CP W_1 (CP K ~ 1 ), and U is a random unitary matrix 
generated according to the Haar measure on U (NK) . It is easy to see that this construction leads to the induced 



measure Pn,k given by (3.6), since the Haar measure on U (NK) is invariant with respect to multiplication by product 
matrices U (N) ® U(K) rep resenting local operations. Thus our results concern not only average entanglement of pure 
states, but in spirit of |32] also averaged entangling power of global unitary operators. 

IV. GENERALIZED ENSEMBLES OF DENSITY MATRICES 

One may obtain a different ensemble of density matrices starting from a real pure state with time reversal symmetry 

K 



<J 



= ^u ik u jk (4.1) 



k=l 



with real Ui k . The only constraint we have now is ^2 ik uf k = 1. Again we may restrict to K > N and since g is 
positive definite we can make the transformation 

u = ^u, (4.2) 
Now the matrix-delta function of the real symmetric matrix transforms according to 

6{y/g{l - uu T )y/g) = detg(- N - 1] / 2 5(1 - uu T ) , (4.3) 

such that 

P(q) oc 9(g) S(Trg - 1) detg (K - N - 1 ^ 2 (4.4) 
and the distribution of eigenvalues is given by 

P« (A lf . . . , \ N ) oc 5(1 - £ A,) J] \\ K ~ N ~ 1 ^ 2 6(Ai) J] \K — Aj I . (4.5) 

i i i<j 

Thus Pn.k — Pn^k corresponds to unitary symmetry with (3 = 2 and P^ 1 -* to orthogonal symmetry with = 1. Both 
are special cases of 

i i i<j 

which can also be given a meaning in the case of symplectic symmetry /3 = 4 (in that case each Kramers degenerate 
eigenvalue is counted once and rescaled). In the general case of arbitrary real > the normalization constant reads 

c (ff) T(KN0/2)(T(l + 0/2)) N 

N ' K Il?Jo n(K - i)/3/2)r(l + (N- j)0/2) ■ 
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It can be derived using the map to the Laguerre ensemble and using Selberg's integral These ensembles of 

density matrices may be constructed by generating real {[3 = 1) or complex (/3 = 2) random Gaussian matrices Ap 
and defining random density matrices by g = ApAU (ti ApAV). In the case of symplectic symmetry ((3 — 4) one has 
complex matrices Ap but then takes only the selfdual part of p. In fact, any density P(A), in the space of complex 
density matrices of size N generates in this way a certain ensemble of density matrices of this size. 

Let us mention that non-negative matrices AA> are sometimes called random Wishart matrices. The joint probabil- 
ity distribution of their eigenvalues (i.e. the distribution of singular values of A) was analyzed in f33| , |34j| . Real random 
matrices A±A[ may describe the matrix of correlations between time series of different stock prices in the presence of 
noise. This fact explains a recent interest in their spectral properties (or in the distribution of singular values of real 
random matrices A\), from the point of view of mathematical finances |p5|-[3^] . Thus the theory of random matrices 
provides an unexpected link between the computation of the mean entanglement of a certain ensemble of composite 
quantum states and the estimation of the financial risk ]38| . 



V. CONCLUDING REMARKS 



We described various ways of defining ensembles in the space of density matrices of finite size N. While it is natural 
to take eigenvectors distributed according to the Haar measure on U(N), there are several possibilities of introducing 
the measure into the (N — 1) dimensional simplex of eigenvalues. Besides measures related to random vectors (of 
orthogonal or unitary matrices), one may define measures related to distances in the space of density matrices (say, 
Hilbert-Schmidt measure or Bures measure). 

An alternative way consists in assuming certain measures in more dimensional spaces (e.g. in the space of pure 
states in NK dimensions or in the N 2 dimensional space of square matrices of size N) , and considering the measures 
induced by projection onto the simplex. We generalized the results of Hall [Q deriving an explicit formula for Pn,k 
and p roved that the measure induced by symmetric partial tracing, Pn,n coincides with the Hilbert-Schmidt measure 
(2.12). We succeeded in computing some averages with respect to this measure including an explicit formula for 
the Shannon entropy of the reduced density matrix, equal to the mean entropy of entanglement, averaged over the 
manifold of pure states of the composite N x N system. 

Moreover, we demonstrated a link between the H-S measure determining the spectra of the reduced density matrices 
and the singular values of the random matrices distributed according to the Ginibre ensemble. As a by-product we 
demonstrated that N 2 random complex (real) numbers drawn with respect to the normal distribution, which constitute 
a non-Hermitian (non-symmetric) random matrix, can be considered as a vector of the unitary (orthogonal) random 
matrix of size iV 2 , drawn according to the Haar measure on U(N 2 ) (on 0(N 2 )). It would be interesting to find an 
ensemble of random matrices corresponding to Bures measure ( 2.15| ), to compute the normalization constants C' N for 
arbitrary N and to evaluate the mean entanglement averaged with respect to this measure. 

It is a pleasure to thank Pawel Horodecki, Marek Kus and Thomas Wellens for inspiring discussions. Financial 
support by Komitet Badah Naukowych in Warsaw under the grant 2P03B-072 19, the Sonderforschungsbcrcich 'Un- 
ordnung und grofie Fluktuationen' der Deutschen Forschungsgemeinschaft and the European Science Foundation is 
gratefully acknowledged. 

Note added: After this wor k was completed we learned form M.J.W. Hall about related papers published on this 
subject: the distribution (3.6) was derived by Lloyd and Pagels p9| , while the formula for the mean entropy averaged 
over the induced measures was conjectured by Page [Q and proved by Sen |4lJ , (see note added to |Q, preprint ArXiv 
quant-ph/9802052). We are grateful to M.J.W. Hall for helpful correspondence. 



APPENDIX A: RESCALING OF RANDOM GAUSSIAN VARIABLES 



We are going to analyze the joint probability distribution of independent random numbers y[, rescaled according to 

y'i -> vt = ■ (Ai) 

We assume that each random number is drawn according to the same probability distribution P(y') defined for 
positive y' . It is easy to see that the joint distribution of the rescaled variables, P(yi, ...,yjv), defined in the interval 
I = [0, 1], is symmetric with respect to exchange of any variables, yi < — ► yk- It is worth to point out that the uniform 
initial distribution, P{y[) = 1 for y' <E [0, 1] does not produce the uniform distribution of the rescaled variables, but 
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P(yi....y N ) = 6(1 - y l )/(N(max{y l }) N ). In the simplest case, N = 2, it gives P(y) = 1/(2(1 - y) 2 ) for y < 1/2, 
(where y = y\ = 1 — y%) and symmetrically for y > 1/2. 

Let us now assume that the random numbers are given by squares of Gaussian random numbers, y' i = x\ and 
P(xi,X2, ---Xn) oc exp(— J2i x i)- We want to calculate the distribution of the rescaled variables yi — x 2 /J2i x i : 

P(y%,y2, — 2/iv) oc dxi...dxpfexp(— 

E^n^-^/E^)- ( A2 ) 

i i i 

Introducing a i5 function for the constraint t = ^\ x\ we have 

P(y u y 2 ,...y N )<x [ die"' ( dx 1 ..Ax N 5(t - E TT S (Vi ~ x iA) ( A3 ) 
Jo J i 

If we now rescale — ► a^Vi we obtain 

P{ V i,V2,...Vn)<x J°° dte-H- 1+N ^ J d Xl ...dx N 6(l - E IJ% " ^) ( A4 ) 



which yields the Dirichlet distribution (2.9) with s = 1/2. 

This result may be generalized for the sum of (3 squared random numbers drawn according to the normal distribution, 
y'i = ^2j=i x i ji ^ e distribution of the sum is given by the xf} distribution with (3 degrees of freedom. In an analogous 
way we obtain that the joint distribution of y^s rescaled as in (|Al[ ) is given by the Dirichlet distribution 

N 

P(yx,V2, ...,y N ) « 5(1 - E V*) II wf /2_1 > ( A5 ) 



with s = (3/2. The most important case, f3 = 2, shows that the distribution of squared components of a rescaled 
complex vector of Gaussian random numbers is uniform in the (N — 1) dimensional simplex. 

Above results find following simple applications. Real (complex) random vectors, i.e. columns (or rows) of orthog- 
onal (unitary) random matrix may be constructed out of independent random Gaussian numbers (real or complex, 



respectively) rescaled according to Eq. (Al). This is consistent with the well known fact that the distribution of 
squared components of eigenvectors of a GOE (GUE) matrix are given by the Xp distribution with /3 = 1 (or 2) 
degrees of freedom |l6|]. Thus certain averages, one wants to compute averaging over the set of pure states (say, 
entr opy, participation ratio, etc), may be replaced by averages over the simplex in R N_1 with the Dirichlet measure 



( |2.9D with s = 1/2 or s = 1. 

Furthermore, a complex square random matrix of size N pertaining to the Ginibre ensemble after rescaling 



ip' — > ip = ip' /y trip'^tp' represents a random pure state drawn according to the natural measure on 

cp n -i_ a similar 

link is established between rescaled rectangular matrices of size K x N and the elements of CP KN ^ 1 . Analogously, 
normalized non-symmetric real random matrices with each element drawn independently according to the normal 
distribution are equivalent to random vectors distributed with respect to the natural measure on the real projective 
space, RP N2 -\ 
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